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Abstract 



Ph ' We investigate the inhomogeneous Galton-Watson processes with 

^S^ . immigration, where pn the offspring means in the n*^ generation tends 

j^ I to 1 . We show that if the second derivatives of the offspring generating 

functions go to rapidly enough, then the asymptotics are the same as 
in the INAR(l) case, treated in [4]. We also determine the limit if this 
assumption does not hold showing the optimality of the conditions. 
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O ; 1 Introduction 

Let Xq = and consider the following inhomogeneous Galton-Watson pro- 
cess with immigration: 

i=i 

where {^n,j,£n '■ n,j G N} are independent nonnegative integer valued ran- 
dom variables such that {C,n,j : J G N} are identically distributed. If the 
offspring distribution is Bernoulli distribution, that is each particle either 
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dies without descendant or leaves exactly one descendant, we obtain the so- 
called first order integer-valued autoregressive (INAR(l)) time series. We 
assume that the process is nearly critical, that is 

E^n,i = Pn t 1 as n — 7- oo. 

We note here that in the followings any non-specified limit relation is meant 
as n — )■ oo. 

The theory of branching processes is used to model the evolution of pop- 
ulations whose members live, reproduce and die independently of each other. 
Its first appearance was motivated by the problem of extinction probability 
of the family names in the British peerage. The problem was addressed by 
Francis Galton in 1873 and solved by Henry Watson in 1874. Since then 
the theory is developing to model more complex systems, and now branching 
processes play an important role in models of genetics, molecular biology, 
physics and computer science. Therefore we do not even try to give a com- 
prehensive bibliography. As a main reference on branching processes we refer 
to the classical book of Athreya and Ney [1] . For some recent application of 
INAR models and branching processes see [4] . 

The aim of the present paper is to investigate the asymptotic properties 
of nearly critical Galton- Watson processes with immigration under general 
offspring distribution, and thus extend the results in the INAR(l) case by 
Gyorfi, Ispany, Pap and Varga [4], by dropping the restrictive condition of 
Bernoulli offsprings. 

Section 2 contains all the results and the discussions. In subsection 2.1 
we investigate the case which is parallel to the results in [4]. Here we assume 
that the variance of the offsprings tends to with a prescribed rate, and 
thus in a clear sense the offsprings are 'almost' Bernoulli random variables. 
The methods here are similar the ones in [4] , however we emphasize that the 
proofs are necessarily more difficult, since we do not have a closed form for 
the generating functions even in the case of Bernoulli immigration. Subsec- 
tion 2.2 deals with a significantly different case, when the second derivatives 
also contribute to the limit distribution, and so the proofs need new ideas. 
Subsection 2.3 contains the case, when p„ — )■ 1 very fast. Finally, in subsec- 
tion 2.4 we consider the case of linear fractional generating functions. This 
example is very important because of the explicit computations. All the 
proofs are placed in section 3. 



2 Results and discussion 

Introduce the generating functions 

Fn{x) = Ex^", Gn{x) = Ex^"-\ Hn{x) = Ex'" , X G [0, 1]. 
Using the branching property we obtain the basic recursion for F„: 



FJx) 



E 



X 



Z^j^—l ^n,i+£rj 



E 



E 



X 



Z^i — l ^n.x 



X. 



n—1 



E [G'„(X)^"-] Hn{x) = Fn-l{Gn{x))Hn{x). 



Introduce the following notation: G„_|_i „(x) = x and if Gj+i_„(x) is defined 
for j < n then 

Gj,n{^) = (GjO . . . O Gn) (x) = Gj{Gj+i^n{x)). 

With this notation the induction above gives the formula 



F„(a;) = n^.(G,+i,n(:r)). 



(1) 



According to the continuity theorem for discrete random variables ([2] 
p. 280) for proving a limit theorem we have to show that Fn{x) converges as 
n — )■ cxo for all x G (0, 1), and the limit function is the generating function 
of the limit distribution. Since for fix j the function i7j(Gj+i „(a;)) ~ 1, we 
introduce the corresponding generating function 



^n(a:) = J]e^^(^^+^-(^))-i = exp \ ^(i7,(G,+i,„(x)) - 1) L 



j=i 



.i=i 



which is easier to handle, because of its exponential form. This is a kind of 
accompanying law of X„. So for proving a limit theorem we have to check 
the following two conditions: 

(a) Fn{x) — Fn{x) — 7- for all x G (0, 1), and 

(b) the limit lim„^oo-^n(2;) exists for all x G (0, 1), and the limit function 
is a generating function. 



Also note the intuitively clear fact, that 

G'i+l,n(l) = Pj+lPi+2 ••■ Pn = P[j,n]- (2) 

This is because (?„„(!) = ^'^(1) = p„ and 

and so by induction we have (2). 

2.1 Poisson and compound Poisson limits 

In this subsection we investigate the case, when the second derivative of the 
offspring generating function goes to so fast that it does not appear in the 
limit distribution. So the results here are the analogs of the ones in the case 
of Bernoulli offsprings in [4] . 

In the followings Poisson(A) stands for a Poisson distribution with pa- 
rameter A if A > 0, and Poisson(O) is the degenerate distribution at 0, while 
Bernoulli(p) is a Bernoulli distribution with parameter p > 0. 

The first theorem deals with the case when the immigration has Bernoulli 
distribution. 

Theorem 1 Let {Xn}nm be a Galton-Watson process with immigration, 
with En ^ Bernoulli (?72„^i). Assume that 

(i) pn < 1, lim„^oo Pn = 1, ErT=i(l ^ Pn) = ^ ^ud hm„^oo Y~^ = ^' 
(ii) lim^^oo Yz^ = A. 

Then 

Xn — > Poisson(A). 

Next we turn to general immigration distributions. Let us denote the 
factorial moments of the immigration distribution by 



mn,k 



E[en{en - 1) ...(£„- A; + 1)] = //«(!). (3) 



(Clearly, this is consistent with the notation Bernoulli(m„,i).) 

The analog of Theorem 2 [4] can be shown: the process has Poisson limit 
even if the immigration distribution is not Bernoulli, only 'close' to it. Since 
the proof is also the same as there, we omit it. 



mn,i 

-71— >-00 1 „ 
^ Pn 


= X, lim„^oo pf = 0. 




Xn — > Poisson(A). 



Theorem 2 Let {Xn}nen be a Galton-Watson process with immigration. 
Assume that 

(i) pn < I, lim„_,oo Pn = 1, Er=i(l - Pn) = oo and lim„^oo jzj^ = 0, 
(a) lim. 
Then 



Before the following results we recall that for a finite measure /x on Z+ = 
{1,2, . . .} the compound Poisson distribution CP(/i) with intensity measure 
H is the distribution which has generating function 

oo ^ 

exp<;^/i{j}(a;-'-l) I, xg[0,1]. 

The analogs of Theorem 4 and 5 in [4] are also true. The proof of Theorem 
3 is the same as of Theorem 4 in [4], and also follows from the stronger 
Theorem 4, so we skip it. 

Theorem 3 Let {Xn}ne'M ^^ '^ Galton-Watson process with immigration. 
Assume that 

(i) Pn < 1; lim„^ooPn = 1, Zl^i(l - Pn) = oo and lim„^oo ^^ = 0' 

(ii) lim„^oo j(7-p„) = \ for i = 1,2,..., J with \j = 0. 

Then 

Xn A CP(/i), 

where p is a finite measure on {1,2,...,J— 1} given by 

1 '^^■'^^ (—lY 

p{j} = - J2 —r^^+^^ j = i,2,...,J-i. 

The following theorem is more general than Theorem 5 in [4] even in the 
case of Bernoulli offspring distributions. And since the proof of Theorem 5 
in [4] is only given by Poisson approximation we give the analytical proof. 



Theorem 4 Let {Xn}nen be a Galton-Watson process with immigration. 
Assume that 



(i) pn < 1; lim„^oo Pn = 1, Zl^i(l - Pn) = oo and lim„ 






0, 



'ii) 


lim """J - 


= ^j for j 


= 1,2,.. 




iii"n-^oo j(i_p„) 


■! < 


ker 


I 






Y, 



where the random variable Y has generating function 



1=1 



l\ 



Note that if A„ = for some n > 2, then A^ = for aU m > n. This 
simple fact follows from the previous two theorems. 

Also notice that the assumption lim sup^^^^ \f\Jn\ < 1 implies that the 
limiting generating function exists for x G (0, 1). This assumption is weaker 
than the one in [4], which is that for all j = 1, 2, . . . the limit 



P{j} 



E 



(-1)' 



A 



•' i=0 



■j+i 



(4) 



exists. However, under this assumption the limit turns out to be a compound 
Poisson random variable with intensity measure /i, that is 



1=1 



l\ 



exp<^^^{/}(x'-l) 



1=1 



This follows easily by Abel's theorem. 

The limit can be compound Poisson even if condition (4) fails. If A^ 
(n — l)!/n, then condition (4) does not hold even for j = 2, but 



1=1 
with 



{x-iy 



x~l 



logfl - u) 



-du 



u 



log(2 - y) 
i-y 



dy = ^p{j}{x^ -1), 



i=i 



,U} = - 



fc=l 



that is the hmit is compound Poisson with intensity measure fi. To see that 
the sequence {A„ = ^"~ }^i can be a hmit in condition (ii) consider the 
foUowing example: p„ = 1 — n~^ and Hn{x) = 1 + {H{x) — l)/n with the 
generating function H{x) = 1 — log(2 — x). 

We do not know whether the hmit in the previous theorem necessarily 
has a compound Poisson distribution. 

2.2 Negative binomial limits 

In the foUowings we investigate the case, when G^(l)/(1 — p„) -/^ 0. In 
contrast to the previous subsection we show that in this case the second 
derivatives do appear in the limit. Here we restrict ourselves to the case 
when the immigration distribution is close to a Bernoulli distribution. 

Note that for condition (a) it is only needed that p„ < Ij Pn — ^ 1 and 
^^i(l — Pn) = oo (see the proof of Theorem 1). Since we always assume 
these conditions we concentrate on condition (b), i.e. -in case of Bernoulli 
immigration- on the existence of the limit 

n 

lim y^mj^i{Gj+i^nix) - 1). (5) 

n— >-oo •'— ' 

We try to compute the derivatives at 1 of the components in the sum 
above. First note that since Gj+j^ „(1) = P[j,n], the first derivative in (5) is 

n n 

and since {(1 — Pj)p[j,n]} form a Toeplitz matrix, we obtain that in order to 
get a limit the asymptotic order of t?2„ i must be 1 — p„. For the second 
derivative we have the recursion 

= G';m^,,M)G',+i,ni^r + g;(g,+i,„(x))g;+i,„(x). (6) 

Substituting x = 1 for j < n (recall that Gn+i,nix) = x) by (2) we obtain 
that 

7 



^11 



and G„„(l) = (^''(1). Induction argument shows that 

n 
^j+l,n(l) = 2^ G'"(l)py^j_i]P[j„], 

SO the second derivative in (5) is 

n n i—l 

j=l i=2 j=l 

As we have seen m„ ~ 1 — p„, and so to get a hmit for the second derivative 
we must have G'^{1) ~ 1 — p„. These heuristic argument kind of shows the 
necessity of the assumptions in the following theorem. We also note that Step 
2 of the proof Theorem 5 shows that the higher derivatives must have the 
same order, namely 1 — Pn- We are not able to calculate these contributions, 
therefore we have to assume condition (iii). 

The reasoning above can be made rigorous to show the following: If 
lim„„^ooG^(l) = a exists and the inhomogeneous Galton- Watson process 
with immigration has a proper limit distribution with finite second moment, 
then necessarily a = 0. This immediately implies that liran^ooGnix) = x, 
i.e. there is no critical nontrivial branching mechanism with finite second 
moment, which can cause a proper limit distribution with finite second mo- 
ment. This result is in complete accordance with Theorem 1 (ii) by Foster 
and Williamson [3] in the homogeneous case. 

We will show that in this setup the limit distribution is the negative 
binomial distribution. A random variable X has negative binomial dis- 
tribution with parameters r > and p G (0,1), denoted by NB(r, p), if 
P{X = k} = {''^L~i^){l -pYp^, A; = 0, 1,2, ..., where the binomial coeffi- 
cient is defined by ( ^^^^ ) = - — ^^^, '" ''"'' . The generating function is 

X I 1-P 



Ex 

1 — px 

The following theorem holds. 

Theorem 5 Let {X„} he a Galton-Watson process with immigration, with 
general offspring and immigration distribution, such that the followings hold: 

(l) pn < 1, lim„^oo Pn = 1, ErT=l(l " Pn) = ^, 



(ii) lim„^oo Y^ = z/ G (0, oo), 
(in) lim„^oo %-r^ = 0, for all s >3, 



(iv) lim„^oo rir = -^ ^^^ lim„^oo fv" = 0- 

Then 

X„ AnB(2A/z/,z//(2 + z/)). 

As we already mentioned, assumption (ii) means that the second deriva- 
tives have a significant role in the limit, while assumption (iii) ensures that 
the third and higher derivatives do not count. The proof of the theorem 
basically lays on determining the asymptotic of G^_,_i „, combined with a rel- 
atively closed formula for the coefficient of f"{g) in -^f{g{x)), see (13). This 
'shows' that to handle the case, when the higher derivatives also count, a dif- 
ferent approach is needed, or at least the calculations become more technical. 
It would be also interesting to extend the results to more general immigration 
distribution, that is to know whether a kind of analog of Theorem 4 remains 
true. 

Finally, we note that if A > is fixed and z/ — )■ 0, then 

lim NB(2A/z/, z//(2 + z/)) -^ Poisson(A). 

It is easy to check that the proof of the theorem remains correct in this case, 
and we obtain that Theorem 5 holds ii u = and in this case the limit 
distribution is Poisson(A), which is exactly the statement of Theorem 2. 

2.3 The case X^^i(l — Pn) < oo 

Let us consider the nearly critical inhomogeneous Galton-Watson process 
with general offspring distribution and general immigration. It was always 
assumed that p„ does not converge too fast to 1, that is X]^i(l ~ Pn) = oo. 
In this section we investigate the case, when this assumption does not hold, 
and assume that 



l[pn = pe{0,l). 



(Note that this is equivalent to the assumption X]^i(l ~ Pn) < oo.) We 
show that in this case the hmit distribution exists under very general as- 
sumptions, but the process does not contain enough randomness as in the 
previous cases: the limit distribution explicitly contains each offspring- and 
immigration distribution. That is p„ tends to 1 so fast, that the process has 
no time to forget the initial distribution. This shows that the right assump- 
tion is indeed Yl'^=ii^ ~ Pn) = oo, as it was investigated in [4]. In fact, the 
case treated in this section is much simpler. 

We state the theorem in the most general setup. 

Theorem 6 Let X„ be the Galton-Watson process with immigration, de- 
scribed above. IfYY^=2Pn = P ^ (0) 1) '^'^'^ Xl^i "^n,i < ^^' then 

X„ —^ Y, 



'-n 



where Y has generating function 

oo 

9{x) = Y[Hj{Gj+i^^{x)), 
i=i 

with 

Gj+i,oo{x) = lim Gj+i.rXx). 

n— >oo 

Let us see an example in the simplest case, when both the offspring dis- 
tribution and the immigration distribution is Bernoulli. In this case we can 
compute the limit generating function. 

Example 1. Let p„ = 1 - ^. Clearly Xl^i(l - Pn) = tt^/G < oo. For 
the product we have 



P[l,n] = P2P3.-.Pn={l-7^][l-:^)---[l--^ 



2y V 32 J ■ ■ ■ V ^ n^ 



(22-l)(32-l)...(n2-l) 



(n\V 



1 ■ 3 ■ 2 ■ 4 ■ . . . ■ fn - 1) • (n + 11 



W)2 



n + l 1 



^ o =P- 



2n 2 

10 



In this case G'j+i^„(a;) = 1 + P[j,n](a^ — 1)) so Gj+i,oo(a;) = 1 + p{x — l)/p[ij]. 
Since Hj{x) = 1 — itlj^i{1 — x), for the hmit generating function we have 



9i^) = n 



1 ^-(I-X) 



n 



1 — ITLji 



J + 1 



:i-a;) 



Choose m^ 1 = (j + l)/j^, so 



9{x) = n 



1 — a; 



Sm TT 



v^ 



a: 



TT 



^ 



X 



The following example is Example 1 in [4], which shows that in a special 
case the limit can be Poisson even in this setup. 

Example 2. Assume that the offspring distribution is Bernoulli(p„) and 
the immigration £„ has Poisson(m„^i) distribution, where X]^i(l ~ Pn) < oo 
and X]^i "^n,i < ^^- ^s before Gj+i,oo(a^) = 1 + p(a: — l)/p[ij], and since 
Hnix) 



gm„,i(a; 1) ^g j^g^^g ^j^g^^ 



g[x) 



exp < (a; 



i)E 



mn,ip 



1 P[l,n] 



and since the sum in the exponent is finite, this is Poisson distribution with 

Eoo rrin IP 

"=1 /'[l.n] 



2.4 The linear fractional case 

The importance of the linear fractional generating functions in branching 
processes is that, that this is basically the only example when explicit com- 
putation can be done. In this subsection we investigate this example in detail, 
which actually helped to find the general form of Theorem 5. 



l-/3s 



The linear fractional generating function has the form /(s) = 1 

where a, /3 G (0, 1), a + /3 < 1. For the first two derivatives we have 



1-/3 



+ 



/'(I) 
/"(I) 



a 



(1-/3)2 
2a/3 
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and these determine the parameters a, /3: 

4/'(l)3 



a 



(3 



(2/'(l) + /"(1))2 
/"(I) 



2/'(l) + /"(!)• 

If /i(s) and /2(s) are both hnear fractional generating functions then so 
is /i(/2('S))) that is if Gn is hnear fractional for aU n, so is Gj^n{s) for aU j, n, 
with parameters 






(2G,„(1) + G,„(1))^ 



/3 - ^^'"^^ 



2g;.ji)+g;;„(i: 



As we have seen G'^.„(l) = P[j-i,„] and ^^^(l) = ELj G'"(l)P[i-i,i-i]P[i,n]- 
The generating function has the form Fn{x) = YYj=i^ji^j+i,n{^))- Assum- 
ing Bernoulli immigration, i.e. Hn{x) = 1 + mn{x — 1), the corresponding 
accompanying generating function is 



Fn{x) = exp<^mj{Gj+i 
= expJ^m, (-- 



Q^j + l,n Oij+l,n-^ 



Pj+l,n ^ Pj-\-l,nX 



Some calculation shows 

Oij+l,n _ P[j,\ 



i-PH.n^ (i+j:..^^,s«3^)(i + j: 






Eventually, we obtained that in the linear fractional case the limit exists 
(under the assumption X]^i(l ~ Pn) = oo) if and only if 

-1 



J™I]^^-^b-.^ 






^^ y^ Gi{^)Pli,n] \ r _^ A G[(l)p[i^ 



,=,+1 2p, y \^ ^^^ 2p, 



exists. It is easy to check that if (ii) in Theorem 5 holds, then so is (iii), that 
is the limit is negative binomial distribution. 

12 



3 Proofs 

We use the continuity theorem for generating functions ([2] p. 280) and Lem- 
mas 5 and 6 in [4] without any further reference. We also frequently use the 
simple facts that 



n--n 



Wk 



k=l 



fc=l 



- X^kfe-^fcl) 



(8) 



fc=i 



for Zk,Wk G [—1,1], k = 1,2, ... ,n, (see for example Lemma 3 in [4]) and 
that |e" - 1 - m| < u^ for |m| < 1/2. 

Proof of Theorem 1. By the assumptions the immigration e„ has 

Bernoulli(m„ i) distribution, that is Hn{x) = 1 + m„^i(a; — 1), in which case 
formula (1) reduces to 



Fn{x) = n [^ + ^^AGj+i^n{x) - 1 



Let us define the generating function 



Fn{x) = JJe"^"-i(^^+^'"(^)-^\ 
j=i 

Estimation (9) gives that if m„,i — > then maxi<j<„mj^i|Gj+i^„(a;) — 1| — )■ 0. 
We want to show that for all x G (0, 1) 

(a) Fn{x) — Fn{x) — )■ 0, as n — > oo, and 

(b) lim„^ooF„(x)=e^(^-i). 

Condition (a) is the easier to handle. Using (8), the inequality |e" — 1 — 
u\ < u"^ and the mean value theorem combined with the monotonicity of the 
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derivative, we obtain 



\FJx)-FJx) 



n 

n 



."^i.iCGj+i.nCa;) 



-i)-n[l + m,-i(G,+i,, 



:x)-l 



i=i 



< J2 e™-^(^^+^-(^)-i) - [1 + m,,i(G',+i,„(a;) - 1)] 



n 



< 



J2^li\Gj+i, 



ix)-l\ 



i=i 



n 

< 5^m|Jx-l|2(G;+i,„(l) 



Using (2), finally we have 

n 
\Fn{x) - Fn{x)\ <\X- M'^^^llpfj,n], 

i=i 

and by Lemma 5 in [4] this goes to under the assumptions of the theorem. 
Now consider condition (b). This convergence is equivalent to 

n 

^mj,i[G'j-+i,„(a;) - 1] ^ A(a; - 1). 
Let X G (0,1). The convexity of the generating functions implies Gn{x) > 

^ — Pn + PnX = 1 + Pn{.X — 1). And SO 

Gn-l,n{x) = Gn-l{Gn{x)) > 1 + p„_i(G'„(x) - 1) > 1 + Pn-lPn{x - 1). 

Induction gives that 

Gj+i^n{x) > 1 + Pj+lPj+2 ■ ■ ■ Pn{x - 1) = 1 + P[j,n]{x - 1), (9) 

which implies the lower bound 

n n 

^mj,i[G'j+i,„(a;) - 1] > ^mj,ip[j,„](x - 1) ^ A(a; - 1). 
i=i i=i 
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For the upper bound note that by the previous estimations G'j+i^„(a;) G 
(1 — P[j,n]) 1) for any x G (0, 1). Again by convexity, for y G (1 — Pij,n], 1) we 
have 

V Plj,n] P[j,n] 

< G.ll-py, „]) + !-- 

P\i,n\ P[j,n] 

= l+^,-„(y-l), 



where 



a _ ^-Gj{l-P[j,n]) _^,,c ^ 



by the mean value theorem. Since p[j^„] — )■ thus ^j,„ — )■ 1 as n — )> oo, 
'&i,n t Pj ■ So we get 

Gj^n[x) = Gj{Gj+i^n{x)) < 1 + ^j,„(G'j+i,„(x) - 1), 

and so induction gives 

Gj+l.rXx) < 1 + '&j+l,n'&j+2,n " ' ' '^n,ni^ - 1) =: 1 + ^b>](x - 1). (10) 

Summarizing we have 

n n 

^^^AGj+iA^) - 1] < ^mj^Aj^n]{x - 1). 
i=i i=i 

Note that in the case of Bernoulh offspring distributions the upper and lower 
bounds are equal. 

So we have to check that under what conditions 

n 

For this, exactly the same way as in Lemma 5 in [4] we only need that the 
sequence bnj '■= (1 — Pj)'^[j,n] form a Toephtz matrix. The only nontrivial 
condition is 

n 

15 



We may write 

n n n 

j=l 1=2 1=3 

+ (1 -p„_l)^„,„ + (1 -p„) 

\Pn '^n,n) ' \Pn—l '^n—l,n)'^n,n ~r ' ' ' 

n n 

+{p2 - ^2,„) n ^'." +pIi ^i'^ 

1=3 1=2 

Note that in the bracket every term is nonnegative, since dj^n = G'A^j^n) < 
G'jil) = pj, because G'j is monotone increasing since Gj is convex. So we 
need that 

n n 

{Pn - '&n,n) + {Pn-l - '&n-l,n)'&n,n H ^ (P2 " ^2,n) Y\_ ^'." + /'l 11 ^''" ~^ ^■ 

1=3 1=2 

By definition we have i!}j^n = G'j{C,j,n), where ^j,„ G (1 — Pij,n], 1), and since 
Pj = G'j{l) the mean value theorem again gives 

p, - ^,,„ = g;.(i) - g;.(4„) = (1 - ^,,n)G';ie,,n) < PMG';ii), 

where ^j„ G (^j>, !)• Therefore we obtain 

n n n (^"(1 \ 

j=i j=i j=i Pj 

which, according to Lemma 5 [4], goes to 0, if ^_ ' — )■ 0. The proof is ready. 



Proof of Theorem 4. By formula (1) for the generating function we have 
to show that for all a; G (0, 1) 

Fn{x) = U^AGj+iA^)) ^ exp 5] ^^TT^^' 
j=l 1 1=1 

Let us define the function 

n 
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Using the estimation (9) and that < 1 — Hj{x) < Tnj^i{l — x), we have 

n 

\F4x)-F4x)\ < J]|e^^(^^+i-(-))-i-//,(G,+i,„(x)) 

n 

< $^(i/,(G,+i,„(x))-l)' 

n 

i=i 

since mj^i/(l — pj) — > Ai imphes m|;^/(l — pj) — )■ 0. So we obtain the 
convergence Fn{x) — Fn{x) — )■ 0, for all x G (0, 1). Therefore what we have 
to show is that 

j=i 1=1 

Let us fix an e > and an a; G (0, 1). There is an /q such that 



tK^x. 



l\ 

1=Iq+1 

By Lemma 6 in [4] we have 



< e and ^^^^yyT^'o^^ ^ '■ 



:ii) 



^^•(^) = E ^(^ - 1)' + ^^vo+i(a:), 



i=0 
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where \R,^i,+i{x)\ < (1 - x)'"+^m,-i„+i/(Zo + 1)!. Thus 



J2 [H,{G,^^,n{^)) - 1] 

n f lo 



E \ E ^(G,^U^) - ly + i?,,o+i(G,+i,n(:r)) 



j=l k 1=1 






1=1 j=l 

lo n 



i=i 



EET {G,^iA^)-iy-Pia^-iy 



1=1 j=i 



E ^7^ E "^Mjm + h + h. 



Using (9) again 






:t(^o + 



(/o + i; 



^«()+l 



by Lemma 5 in [4], thus (11), the choice of /q, shows that |/i| <2e for n large 
enough. In order to estimate I2 we use the inequahties (9) and (10) we have 

(G,+i,„(x)-l)'-p|^.„](a;-l)' 

< I {Gj+i,n{x) - 1) - p[j,n]ix - l)\lp\jl^il - x) 
</(l-a;)'(py,„]-%„])pg], 



i-i 



therefore 

{i-xY 

^ (/-I)! ^ 

1=1 ^ > j=i 

which goes to 0, due to Lemma 5 in [4] and to our assumptions. 



1^2! < E (/_!)! Y'^o.MoM - ^b>])p[,-^], 



Finally for the difference 



n oo / -. w 



i=i 



1=1 



n 



< 



E 

1=1 



[X 



l\ 






E 

l=la+l 



\X 



l\ 



-\l 



IAI + I/2 



where the first and fourth term converge to 0, while the sum of second and 
the third is less than 3 e. Since e was arbitrary the proof is ready. ■ 

Proof of Theorem 5. We separate three cases. First we assume that the 
offspring generating functions are second degree polynomials, that is each 
particle has at most two offsprings. Then we extend the proof when for the 
higher derivatives only assumption (iii) holds. In this two cases we assume 
Bernoulli immigration. In the third case we extend the result for general 
immigration. Step 1 shows the main idea without the technical difficulty to 
handle the higher derivatives. 
Step 1. Let us assume that 



deg G„ = 2 for all n G N. 



(12) 



First we obtain a recursion like (7). What we need is a general form of a 
derivative of a composite function. Let 



hk{x) 



dx^ 



f{9{^))- 



According to Lemma 5.6 in [5] the general formula for the derivatives of a 
composite function is 



1 



da;*^ 



/(,(.)) = A;! 5^ /(^)(,(x)) E HttU 



{ui,...,Uk)m=l 



1 



iMi 



X 



(13) 



where the summation is carried out over all nonnegative integer solutions of 
the equation system: 

z/i + 2z/2 H h kuk = k 

jyi + jy2-\ huk = s. 
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Clearly the coefficient of f'{g) in h^ is g^^\ and induction argument shows 
that 



fe-i 

I y 

ClI I 



E»ii {^^9^^9^^~"\ A; odd, 

coefficient of f"(g) in hk = \ ^ 

Et/©^«^(^-^) + iay(^(t))^ A; even. 

For simplicity introduce 

Wfe,i — \ l( k \ ■f-_k 
2U/2J' ^^ ^ ~ 2- 

Apply this result to Gj^n{x) = Gj{Gj^i^n{x)) for j < n, and note that the 
third and higher derivatives vanishes. Substituting x = 1 we have 

fc/2 

G^lii) = p.Gf^\,ji) + g;'(i) j2 «^.Gj+i,n(i)^f+;;i(i). (14) 

From this, easy induction argument gives the following: 

k/2 n 

Gf!i,nil) = Y.^k, Y. PW-i]G[\l)Gfl,Jl)4t^l{l). (15) 

Until now everything hold in general. Now we use our assumptions to 
prove 

Gfi,„(l) = (^^'"' + «(1)) Py,n](l - PU^n])'-' + Pb-,n]0(l), (16) 

where o(l) — ?■ as j, n — )■ oo. The proof goes by induction, uses the recursion 
(15) and the identity 

(fc + 1)! _ 1 ^^^^' i\ {k + l-i)\ 

Formula (16) is true for k = 1,2, and let us assume that it is true until some 
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k > 2. Then using the induction hypothesis 



„(fc+i) / ^ 



fe+l 
2 



2_^ak+i,i 2_^ P[i,/-i]G'z'(l)G'z+i,„(l)G';_^i,„ (1' 

k + l 
2 " 



i=l 



X 



X 



fc+i 

2 



t\ 






2k- 



G''(l] 



J^ afc+l,iP[i,n] 5^ q'_ x (l-pOP[/,n] 



{k + iy. 



X 



Z/ + o(l) p[j-„](l - py,„]) + py,n]o(l), 



where we used the estimation 

n 

Y^ (i-p/)pm(i-pm) 
l=j+l 



[^ - P[j,r. 



\k+l 



k + l 



+ 0(max(l-p;)p[z,„]). 

J<l<n 



Note that the left side is a Riemann approximation of the integral J{l—y)''dy, 
corresponding to the partition {p[z,n]}r=j+i- So (16) is proved. 
From (16) using also assumption (iv) in the theorem we have 



^.,lG,+i,„(l) = }_^ Y^ [ ^^'^ 

7 = 1 7=1 ^-^ ^ 



+ o(l) Pb>](l-Pi)(l-Pb>]) 



fc-1 



T^V^^ " Pj)Pb-,nF(l) 



i=i 



(* - D'-A.'-' 



)fc-l 



<A;-1)!A 



fe-i 
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So for the sum (5) we have 






j=i j=i k=i 

k=l 

= — -log(^l- (x-l)- 

and Lemma 6 in [4] makes the calculation rigorous. Note that since Gk,n{x) 
is a polynomial the infinite sum above is in fact finite. 
We obtained that the limit generating function is 

exp (-—log fl - (x- 1)^)1 = fl - (x- 1)''^"^ _ / 2T. 



which is the generating function of a negative binomial distribution with 
parameter r = 2\/v and p = z//(2 + z/), as we stated. 
Step 2. Let us weaken the condition (12) on the generating functions. 
From formula (13) induction argument shows that 

s=2l=j + l (,yi,...,,y^)k=l "*■ \ 

We claim that under assumption (iii) of the theorem (16) holds. By Step 
1, this is true for A; = 1,2. Assume that the statement is true until some 
k > 2. The previous case shows that we get the asymptotic from the term 
corresponding to s = 2, and we show that the terms corresponding to s > 3 
are o(p[j„]). Using the induction hypothesis for the term corresponding to 
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s > 3 above we have 



7 11 / \ I 1 "l" \ 

/=7+l i^i,...,;/fe) A;=l \ 



Z=j+1 (;/i,...,i/fe) A:=l 



d^^di 



const py,„] J2 J7Y—;)('^-P'^Pv 



i£jt, Pii^ - Pi) 



s-l 

n\ 



Pb>]0(l)- 



We proved (16), and the statement follows as in the previous case. 

Step 3. Finally, in case of general immigration distribution, according to 

the already proved part what we have to show is 



n^, (Gj+iA^)) - n [1 + ^^M (G.+i,"(^) - 1)] 



i=i 



i=i 



< E 1^^- i^HiA^)) - [1 + ^.-,1 {Gj+iA^) - 1)] 



j=i 



<E^(^.w-)-i)^ 



i=i 



< 



{x-lf 



E "^i.2P[j>] ^ 0, 



where we used Lemma 5 and 6 in [4], and the assumption m„^2/(l — Pn) — ^ 0. 
The proof is complete. ■ 

Proof of Theorem 6. The generating function of the n^^ generation is 



FJx) 



l[H,{G,^^,4x)). 



i=i 



First we show the existence of Gj ^oo{x). Let us fix j, and investigate 
Gj^n{x) as n — )■ oo. Using the definition and the monotonicity of Gj^n and 
that Gn{x) > X for all n, we have 

Gj^n+l{x) = Gj^niGn+l{x)) > Gj^n{x)- 
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Since Gj_„(x) < 1, the limit 

Gj,o:,{x) = lim Gj,„(x) 

exists. Moreover Gj+i^n{x) > 1 — P[j,n](l — x) implies 

G,+i,oo(:r) > 1 - P^}—^. 

Pihj] 

Therefore we can define the function (7 in the theorem. 

Using the inequality above, the convergence rrij^i — )■ 0, the estimation 
Hj{x) > 1 — mj^i{l — x) and that for x small enough 1 — x > e~(^+^)^ we 
obtain 

/9( 1 — X) 

Hj{Gj+i^ooix)) > 1 - rrij^i (l - Gj+i^ocix)) > 1 - ruj^i 

P[id] 



I P[i,j] J 



This easily implies that the function g is continuous at 1, g{l) = 1. 

Next we show that Fn{x) — )■ g{x), for all x G [0, 1]. Introduce the notation 

n 

9n{x) = YlHj{Gj^i^ooix)). 
i=i 
Clearly gn{x) — t- g{x), so we only have to show the convergence gn{x) — 
Fn{x) -> 0. We have 

n 
\gnix) - Fnix)\ < ^ \Hj{Gj+i^ocix)) - Hj{Gj+i.aix))\ 

i=i 

n 
- X]"^i.l \Gj+l,ni^) - Gj+i^ooix)\ 

i=i 

no 00 

j=l i=no+l 

where the first term goes to for every fixed no, while the second one can be 
arbitrary small by choosing no large enough. ■ 
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